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r— s We propose a scheme to realize coherent control of spin states of molecule magnet, Ni4. We 

fvj , introduce transverse magnetic fields with special frequencies. When the frequencies of transverse 

magnetic fields match in some conditions, we obtain dark states in Ni4 molecules. Through adjusting 
the magnitude of magnetic fields, we can obtain any arbitrary superposition of the two ground spin 

^^ , states of Ni4 molecules. 
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The dark state[l| is a specific coherent superposition of ground states that is immune to further 

C/3 ' excitation due to destructive quantum interference. It has played an important role in studies of 

a number of coherent phenomena in physics such as electromagnetically induced transparency 13, 

coherent population trapping [3J, lasing without inversion |j|, slow propagation of light in a medium[5j, 

optical information storage|3, and so on. 

Coherent population transfer^ is another application of dark states. When the dark states condition 

jrt \ is met, via coherent population transfer, one can produce excitation between states of the same parity, 

for which single-photon transitions are forbidden for electricdipole radiation, or between magnetic 

sublevels. The coherent control of quantum states is a central issue in the emerging fields of spintronics 

'■^ . and quantum information processing. Coherent population transfer is an effective method to realize 

Ch ' the coherent control of quantum states. 

O \ Recently, molecular magnets have received much attention because of observed quantum features 

at the mesoscopic level and their potential uses in magnetic storage and quantum computing. Leuen- 

berger and Loss have proposed a scheme to realize Grover's algorithms in molecular magnets such 

as Fes and Mni2jg|. Hou et al. have suggested to perform CNOT operation in dimer of molecular 

^ ' magnets, [Mn4]2.f3| 

In this paper, we propose a scheme to realize dark states in molecular magnet Nizi. pil IllL lia| . so 
that we can obtain an arbitrary superposition of the two spin states of Ni4 molecule. Therefore the 
■.— 1- , coherent control of spin states of molecule magnets is realized. [Ni(hmp)(t-BuEtOII)Cl]4, refereed 

C^ • to as Ni4, consists of four Ni^^(spin 1) magnetic ions and oxygen atoms at alternating corners of a 

^^^ \ distorted cube, with S4 site symmetry. Ferromagnetic exchange interactions between the Ni^^ ions 

^D ■ lead to an 5 = 4 ground state at low temperature, 

■y ' The corresponding Hamiltonian of Ni4 molecule is given by 
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H = H„ + Hi. (1) 

Here Hq is the Hamiltonian without the external field. It can be written as 

Ho = -DSl (2) 

where D is the axial anisotropy constant and Sz is the z component of spin operator. Hi is the 
Hamiltonian due to the interactions with the applied magnetic field 

id 

5^ ; Hi= -.g/iijB • S, (3) 

where /is is the Bohr magneton, g is the electronic g-factor, S is the spin operator, and B is the 
applied magnetic field. 
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Each Ni4 molecule can be modelled as a 'giant spin' of 5 = 4 with Ising-like anisotropy. Every 
eigenstate of Ni4 can be labelled by \m) with rr?, = 4, 3, • • • , —4. With the projection operator \m) (to|, 
the Hamiltonian without the applied external field can be written as 
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where 



= -m'^D/h. 
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FIG. 1: The diagram of energy levels of Ni4 molecules. Here uj = E/Ti represents the energy of spin states 
and Vi {i — ±1, it2, ±3, ±4) are the frequencies of the applied transverse magnetic fields, (t"*" and a~ denote 
the left and right circularly polarized photons respectively. A is detuning. 

To realize coherent control of spin states of Ni4, we introduce time-dependent transverse magnetic 
fields Vmit) = Bm[cos{iyrnt + (j)m)ex ~ sm{iymt + (f)m)ey] and T4i(t) = Bm[cos{iyrnt + 4>m)ex + sin{iyrnt + 
(j}fh)^y], where m = 1, • • • , 4 and ffi = —4, • • • , —1; 0„i and (jjm are phases of the transverse magnetic 
fields; e^ and e^ are respectively the unit vectors pointing along x and y axes. Thus we obtain the 
Hamiltonian due to the interactions with the transverse magnetic fields as 



Hi{t) = -g^J.B < ^ Bm[cOs{iymt + (j)m)Sx - Sin(t'„i + (t)m)Sy\ + ^ B^[coii[v^t + cj),n)Sx + Sin(l/,ftt - 

l^ 771^1 m^ — 4 

l^m— 1 ni——4 ) 

The transverse magnetic fields described by the first part in Eq.Q rotate clockwise and produce left 
circularly polarized a~ photons. Absorption (emission) of a~ photons gives rise to Am = — l(Ar7i — 
+1) transitions of spin states. While the transverse magnetic fields described by the second part in 
Eq.Q rotate anticlockwise and produce right circularly polarized a^ photons. Absorption (emission) 
of (7+ photons gives rise to Arn = -|-1(Ato = —1) transitions of spin states. 
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There is a key point needed to mention. In order to eliminate the excitation states in the process, 
we choose the special frequencies of the transverse magnetic fields with appropriate detuning from 
the gaps between spin states. Here we choose the frequencies i'i{i = — 4, — 3, — 2, — 1, 1, 2,3,4) of 
the transverse magnetic fields far-detuning from energy gaps of spin states, such as U3 — Wi — v^ — 
3A, a;2 — (1^3 — i^3 — —A,lui — uj2 — i'2 — — A, wq — 1^1 — i^i — — A, wq — t^'-i — i^-i = —A,tu-i—uj-2 — i^-3 = 
—A, uj-2 — LO-3 — v-i — —A, w_3 — LO-i — v-i — 3A and Vi — v-i {i — 1, 2, 3, 4), where A is a finite 
detuning. The detail of the choosing of the frequencies Vi is shown in detail in Figure Q] In Reference 
[l^. measured EPR spectra on the Ni4 molecules reveal that values of D are in the range of 0.72 — 1.03 
K. So, the energy gaps between different spin states are ^ 10~^'^J. We estimate that the frequencies 
of the transverse magnetic fields are in the range of 100 '--^TOOGHz. 

Next we calculate the quantum amplitudes for the transitions induced by the time-dependent trans- 
verse magnetic fields by evaluating the 5-niatrix perturbativelyjg- The 5'-matrix can be expanded 
in the perturbation series in power of the interaction Hamiltonian H[(t) = e^^°*'/^Hi{t)e~^^°''-/^ as 
S = X]^o '5'''"''- The nth-order term of the perturbation series of 5-matrix is expressed by 

/ „■ \ " r+00 j-ti j-t„-i 

S^"'^ = -- / dh dt2--- dt„Hiih)H(it2) ■ ■ ■ Hiitn). (6) 

V '^/ J —oc J — 00 J —oc 

In the Hq representation, the ^'-matrix is expressed by Spa = X^^o '^gJ ^ which means the transition 
amplitude from state a to state /3. After a straightforward calculation, we find that the perturbation 
terms below the fourth order vanish. We also obtain the nonzero fourth order perturbation terms as 
follow, 

(4) _ ie-^i^^+'P^+'Ps+Mnl^n, /•- 

,-p-i(0-i+0-2+0-3+0-4) rr^"'" O- r°° 
SHU = MA3 "'"" " LJ'- («) 

where fli — g^sBi/h [i ~ — 4, — 3, — 2, — 1, 1, 2, 3,4). The higher order perturbation terms are neg- 
ligible compared with S'q j and S^ I4, so we shall not consider the contribution of the higher order 
perturbation terms in our approximation. 

We can interpret the above perturbation series of S'-matrix as the below explanation. Because 
we choose the appropriate frequencies of the applied transverse magnetic fields with detuning, the 
transition rates between some spin states are so small that they can be neglected. Only the transitions 
between spin states |4) and |0) and between spin states | — 4) and |0) are obvious since the frequencies 
of the applied transverse magnetic fields match the gaps of these spin states in four-photon process 
as shown in Figure ^ These transitions are four order processes, i.e. when the transitions happen, 
four photons are absorbed or emitted. If the molecules only occupy the spin states |4) and | — 4) 
initially, so, after the applied transverse magnetic field are added, only the populations of the spin 
states |4), | — 4) and |0) are nonzero. Therefore, Ni4 molecule can be considered as a three-level atom. 
Suppose that we have an effective Hamiltonian -fff^f in the interaction picture about the three- level 
atom model and that the S'-matrix of it is equal to that of Hi. Thus, with this approximation, we 
obtain the effective Hamiltonian in the interaction picture as. 
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i?i^^ = -^[f7,(e-^^»|G)(4|+e^*"|4)(0|) + f7fc(e-'*''|0)(-4|+e*|-4)(0|)] (9) 



where (pa — X)m=i 't'm and 0b — X]m=-4 't'm] ^a and Jlf, are the effective Rabi frequencies written as 



48A3 



^'^ = %^' (10) 
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We consider the Ni4 molecule as a three-level atom described by the effective Hamiltonian 0. For 
convenience, here we set </);, = and 4>a — (t>b ^ ±7r. It is straightforward to verify that the following 



combinations of spin states |4), | — 4), and |0) are eigenstatcs of this Hamiltonian [tIIT^: 



|vI/+) = ^(|O)-sin0|4)+cos0|-4)), 

l^o) =cos6l|4) +sin6l| -4), 

|vI'_) = i=(|o)+sin0|4)-cos0|-4)), 



(12) 
(13) 
(14) 



where cos 6' = flij/ ^yTf^ 
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sin 9 = fla/ y^^'^ + ^l; is mixing angle. The corresponding eigenvalues 
0. Here, the wave function |^o) is the so-called dark state. From 
Ea. (|13|l . we can see that in the dark state the populations of the spin states |4) and | — 4) do not vary 
with time and no excitation to the spin state |0) happen. The population is trapped in the lower spin 
states. Although the transverse magnetic fields exist, there is no absorbtion due to the destructive 
quantum interference between the two transitions. 




90° 

CD 45° h 



^ 0.5 



(a) X \ X \ 




(b) /^\ 




(c) |4> \ / '""^^ 



time 



FIG. 2; The diagram of coherent control of spin states |4) 
angle; (c) the population of spin states 
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When the Rabi frequency Of, is nonzero and fla — 0, i.e. the mixing angle 9 — 0, the spin state |4) 
is right the dark state |\E'o)- In this condition, the gaps between the dark states and the eigenstates 
|^+) and |vl/-) are ±|nf,, which are nonzero. We choose the spin state |4) as the initial state. At 
first we adiabatically increase Qb from zero to a finite value and keep Qa zero. Then we change the 
Rabi frequencies Oa and fib as we want via varying the amplitudes of the transverse magnetic fields, 
so the mixing angle 9 will vary continually with time. Therefore the state can be written as 



l*o(i)) 



nb{t)\4) + na{t)\ - 4) 
^nait) + nb{t) 



(15) 



where flait) and flb(t) are time-dependent functions that vary adiabatically with time. We note that 
any arbitrary superposition of the two ground spin states |4) and |— 4) can be prepared by appropriately 



tailoring the shapes of pulse of the transverse magnetic fields. That is, we can coherently control the 
spin states of molecular magnets by adiabatically changing the Rabi frequencies. 

In order to clearly explain the process of coherent control of spin states, we take two special examples 
that are showed in Figures Inland 13 respectively. In the case showed in FigureEl the initial state is |4). 
In the beginning, the Rabi frequency flf, adiabatically increases from zero and fla keeps zero. It easy 
to verify that the initial state is right the dark state at the moment. Because only the population 
of the spin state |4) is nonzero, the mixing angle 9 is zero. After fib increases to some value, we 
can decrease ilf, and increase O^ from zero adiabatically . When Qb vanishes, f2a arrives to a finite 
value. Simultaneously, the mixing angle 9 changes from to 7r/2, so the coherent population transfer 
is realized from |4) to | — 4), i.e., one can coherently control the evolution of the spin state from |4) 
to I — 4) by this scheme. In Figure 13 we show another case that the coherent and adiabatic evolution 
of the spin state from |4) to -75(|4) + | — 4)) is realized. In this case, firstly, we increase ilb to a finite 
value and fla is zero as the previous case. Next, we keep fib invariant and adiabatically increase Qa 
to the same value from zero. Finally, we decrease ila and r2& to zero simultaneously. The last two 
steps are different from that of the previous case and the mixing angle varies from to 7r/4. In the 
end, the populations of the spin states |4) and | — 4) are both 1/2. From FigureOwe can see that the 
coherent control of spin states from |4) to -^^(K) + | — 4)) is realized by our scheme. 
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FIG. 3: The diagram of coherent control of spin states |4) 
mixing angle; (c) the population of spin states 
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(|4) + I — 4)). (a)the Rabi frequencies; (b) the 



Until now, we have not yet discussed the adiabatic condition for evolution of the dark state in this 
paper. In order to interpret adiabatic evolution of the dark state, we take the case shown in Figure 
121 as example. In Figure 01 we draw the variation of eigenvalues of the three eigenstates with time 
in the condition shown in Figure (21 At the initial time, all the transverse magnetic fields are not 
added, |4) is just the dark state |^o) and the three eigenstates are degenerate. In interval I, Q-b is 
added and fia keeps zero, the energy gaps between the dark state and the other eigenstates appear 
but the eigenvalue of the dark state |\l/o) remains unchanged. During interval II, both l^a and Ob are 
nonzero, the dark state has |4) and | — 4) components and the energy splitting of eigenvalues E±. and 



6 

Eq is largest. In interval III, the Rabi frequency Qb is zero and fla gradually decreases to zero, the 
energy splitting of eigenenergies gradually disappears. From above we can see that there always exist 
gaps between eigenvalues E± and Eq in our operation time. In other cases besides the example, the 
processes are similar. Thus, if we choose the spin state |4), which is right the dark state in interval 
I, as the initial state, it is possible to evolve adiabatically the dark state to our requested quantum 
state by adding special transverse magnetic field pulses. 

In order to realize adiabatic evolution of the dark state, it is necessary that the amplitudes of 
transverse magnetic fields are large enough. If the coupling is insufficient, i.e. the Rabi frequencies 
are too small, the practical quantum state can not follow the evolution of the dark state |\E'o) and loss 
of population due to nonadiabatic transfer to the states |vl'^)and |^-) may occur. The Hamiltonian 
matrix element for nonadiabatic coupling between the state |*o) and the states |^+) or |^_) is given 
by (*±|*o) iij. Thus, the adiabatic condition is 



;*±i*c 



< 



\E±-Eo\ 



(16) 



When the transverse magnetic field pulses have a smooth shape, a convenient adiabaticity criterion 
may be derived from Eq. Hl()|) by taking a time average of the left-hand side, (^±|^o)aii = tt/2At, 
where At is the period during which the pulses overlap. This average value should not exceed the 
right-side, \E± — Eq\ /h — \/^l + ^If^- ^he adiabatic condition is p] 



^av^T > 10 



(17) 



where flav is the average Rabi frequency of Qa and fib, and the number 10 is obtained from experience 
and numerical simulation studies. 

The experimental work |llj | shows that, in Ni4 molecules, the longitudinal relaxation time is long, 
and its order is ~ s. Here we assume that At has the same order with the relaxation time. From this 
relaxation time and Ea. (|17|) . we obtain the restriction for the average Rabi frequencies is Qav > 10s~^. 
In our scheme, it is proper that we choose the detuning A as IGHz. Thus, from Ea. pO|l or l|ll|l . we 
know that J7i > 2.6 x 10^s~^, where Cli = giisBi/h is the average value of fli = g^BBi/h. Here we 
introduce the average transverse magnetic field Bi. Reference [l3| shows g ~ 2.2. Finally we obtain 
the amplitude of the transverse magnetic fields required in our scheme as Bi > 8.4 x 10~''T. In our 
scheme, the precise matching between the sum of frequencies of the left circularly polarized transverse 
magnetic fields X]i=i ^i ^^'^ ^^^ sum of frequencies of the right circularly polarized transverse magnetic 
fields J2i=-4 ^i is required, but a finite detuning between X^j^j^ Vi and the gap of the spin states |0) 

and |4), wq — ^^4, or between X)i-4 ^ ^"^^ ^^^ S^P o^ ^^^ spin states |0) and | — 4), wq ~ ^^-4, does 
not prevent the dark resonance. From the above estimation, we can see that the conditions required 
in our scheme are accessible experimentally. 
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FIG. 4: The variation of eigenvalues of the three eigenstates of the Hamiltonian with time in the condition 
presented in Figure |5| 



In conclusion, using the method suggested by us, one can prepare any arbitrary superposition of 
the spin states |4) and | — 4) from the ground spin state |4) or | — 4). In molecular magnets, there is 
a potential barrier between the spin states with positive magnetic numbers and with negative ones, 
as is shown in Figure ^ In the absence of longitudinal magnetic field, each molecule has a double, 
Kramers degeneracy in its ground states. In Ni4 molecules, they are |4) and | — 4). The transition 
from one ground state to the other is possible only if the molecule jumps over the potential barrier or 
tunnels through it. Moreover, it is difficult to coherently control the transition caused by jumping over 
or tunnelling through the potential barrier and these schemes are not spin-state-selective. However, 
our scheme is spin-state-selective and the process is coherent. Therefore, Through adjustment of the 
magnitudes of the transverse magnetic fields, any arbitrary superposition of the two spin states of 
Ni4 molecule can be obtained, i.e., the coherent control of spin states in molecular magnet, Ni4, is 
realized. 
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